Self-consistent local GW method: Application to 3d and Ad metals 
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Using the recently developed version of the GW method employing the one-site approximation and 
self-consistent quasiparticle basis set we calculated the electronic structure of 3d and 4d transition 
metals at experimental atomic volumes. The results are compared with traditional local density 
approach and with experimental XPS and BIS spectra. Our results indicate that this technique can 
be used as a practical starting point for more sophisticated many-body studies of realistic electronic 
structures. 



I. INTRODUCTION 

The electronic properties of many solids are described 
reasonably well in the local density approximation (LDA) 
of the density functional theory (DFT)j* which is es- 
sentially a single-particle theory representing a cardi- 
nal simplification of the original many-body problem. 
However, the LDA becomes inadequate in materials ex- 
hibiting strong spatial correlations between different elec- 
trons in atoms with localized orbitals, temporal correla- 
tions (retardation) associated with Coulomb screening, 
etc. The corresponding difficulties are extremely hard to 
overcome within the DFT without sacrificing the first- 
principles approach. 

One of the methods used to correct the deficiencies 
of LDA in the description of excited states in metals and 
semiconductors is the GW approximation (GWA)£A The 
GWA usually improves the band gaps in semiconductors 
and insulators^ in metals it may provide information on 
the quasiparticle lifetimes and renormalization which is 
absent in the DFTi&4 However, in all treatments used 
until recently, the GWA was employed in a non-self- 
consistent fashion, by using unrenormalized Green's func- 
tions based on the Kohn-Sham wave functions obtained 
in LDA. Such an approach is easier to implement, but 
it is internally inconsistent, does not correspond to any 
Luttinger-Ward functional 5 and hence violates basic con- 
servation lawspi and produces different results depending 
on the approximation used to solve the Kohn-Sham equa- 
tions. On the other hand, studies of the homogeneous 
electron gas have shown that GW self-consistency wors- 
ens the agreement with experiment at typical metallic 
densities;^ highlighting the limitations of GWA which is 
strictly correct only in the high-density limit. 

Recently, two self-consistent realizations of the Green's 
function method were demonstrated. One of them£ was 
tailored for transition metals and employed the one-site 
approximation for the self-energy which is based on the 
assumption that dynamical screening in d- metals is well 
described at the intraatomic level. Some spin-selective 
diagrams beyond the GW set were also included. The 
results for Fe and Ni were quite reasonable, although no 
improvement was obtained compared to LDA. While the 
choice of GW diagrams is unjustified for the homoge- 
neous electron gas at typical metallic densities and like- 
wise for semiconductors (the diagrams that are left out 



do not contain any small parameter), it was noted that 
the situation may be better in transition metals, because 
high orbital degeneracy provides an additional enhance- 
ment to the diagrams with the largest number of closed 
electron loops favoring the GW set £ In addition, for tran- 
sition metals the "local" (one-site) approximation for the 
self-energy is quite reasonable due to the short screening 
length and greatly facilitates self-consistent calculations. 

Another realization 9 using GW set with the full In- 
dependence of the self-energy was applied to elemental 
semiconductors. It was found that self-consistency and 
accurate treatment of core electrons improve the agree- 
ment with experiment for the band gaps in Si and Ge. 

While the adequacy of GWA for the studies of metals 
and semiconductors has not been firmly established from 
the point of view of the many-body theory, this method 
may be regarded as one of the possible steps toward a 
consistent Green's function-based scheme. Therefore, it 
is important to establish the accuracy of this approxima- 
tion for different materials, and in particular for transi- 
tion metals where, as noted above, there are reasons for 
it to work better than in the homogeneous electron gas. 
To this end, in this paper we calculate the conduction- 
band spectra for all elemental 3d and Ad metals using the 
self-consistent GW approach of Ref. and compare them 
with the observed ones and with those given by the LDA. 



II. SELF-CONSISTENT GW METHOD 

The technique used in this paper was introduced in 
Ref. 0- Here we describe some points in more detail. 

The self-energy S = S&/8G is obtained from the 
Luttinger-Ward generating functional $ defined by the 
set of skeleton graphs^ The Hartree diagram gives the 
local contribution Vk(r), and the exchange diagram con- 
tributes E x = V(r — r') J ImG(r, r', e)de/Tr where V is 
the Coulomb potential. The total contribution of the 
remaining GW sequence (the "correlation term" ) is 

£ c (r,r',e) = - J G(r, r', e - uo)V{v - n) 

x n(ri,r2,cj)W(r 2 y,w)dri<2r2 — (1) 

2m 

Here W is the effective (screened Coulomb) interaction 
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defined by the integral Dyson equation 



W = V- J V(r-r 1 )Il(r u v 2 ,e)W(r 2 ,r , ,e)dr 1 dr 2 (2) 



and the polarization operator 



n(r, r', w) = J G(r, r', e)G(r', r, e + w)de. 



(3) 



Contour of integration in Eqs. CJ,© is directed along 
the imaginary axis and embraces the cut on the real axis 
from the Fermi energy Ep to the external energy. The 
GW approximation may be modified by the inclusion of 
vertex corrections in the polarization operator © which 
may improve the results in some cases. 

The calculations are drastically simplified by the use 
of the one-site approximation 8 (OS A), in which the self- 
energy is calculated only on one lattice site neglecting all 
non-diagonal maxrix elements connecting different lattice 
sites. Thus, in OSA the self-energy depends on energy 
and on the coordinates r, r' belonging to the same unit 
cell. In order to implement OSA we first need to choose 
an appropriate on-site basis. Inspired by the good de- 
scription of the band structure of densely packed solids 
achieved in the atomic sphere approximation within the 
linear muffin-tin orbital method (LMTO-ASA), we used 
a very similar, minimalistic approach. For the basis set 
we also use just one function per each angular momen- 
tum I and its projection m (below we denote L = Im), 
along with its energy derivative. The radial basis func- 
tions 4>i(r) satisfy the "Schrodinger-type" equation 



A, 



-Vh-Si(q) 



Mr) = 



(4) 



where A; is the radial part of the Laplasian, and is 
an integral radial operator whose kernel is obtained from 
Re £ by projecting onto the I subspace: 

Mr,r',e) = ^ 1 Yl f Y L (r)TteX xc (r,r' ,e)Y L (r')dodo' 

ra 

(5) 

where we denoted Yi xc = T, x + E c , Yl are the spherical 
harmonics, and integration is over the directions of r,r'. 
The operator £j may be represented as the sum of a 
local part nf\r) similar to an external potential and 

a nonlocal part £ ; which (after acting on 4>i) gives a 
linear combination of <$>y with V ^ I. 

Similarly to the LMTO-ASA method, the solutions of 
Eq. Q are only found at one energy e/ for each I which is 
fixed at the center of gravity of the corresponding band. 
We may safely discard the imaginary part of the self- 
energy operator because it is small in the vicinity of Ep 
where t\ is usually chosen. Thus, our radial basis func- 
tions are real. The non-local equation is solved by 
iterations. If <f>f is the solution after n iterations, then 
the solution on iteration n + 1 is obtained using auxiliary 



functions /; and gi defined as 



Q + ^-%-S{°(r) 



Mr) = t^ie^f, (6) 
9i(r) = (7) 



according to 0" +1 (r) = //(r) + Agi(r), where the coef- 
ficient A is found by normalizing 0™ + . Just as in the 
LMTO method, in addition to 4>i we compute its energy 
derivative (f>i = dipi/de: 



A; 



k + 



(8) 



To stabilize the solution of the Schrodinger equation 
|@J|, w e added the exchange with the nearest-neighbor 
cells to Yi xc . The corresponding matrix element was sub- 
tracted from Green's function © below. 

With the on-site self-energy H xc (r, r', e) defined for r, r' 
within the single unit cell we can calculate its matrix 
elements between the muffin-tin "eigenfunctions" Xk.v, 
which in turn are the linear combinations of 4>i and 4>i- 
Due to the k-dependence of Xkv the self-energy also ac- 
quires k-dependence, just as the local potential in LDA. 
Finally, on-site Green's function is found by integration 
using the tetrahedron method: 



C(r r> f) - V ^)xtt(r' 



(9) 



where Xku anc ^ xti> are right and left eigenvectors, and 
Aki, the eigenvalues of the non-Hcrmitian operator H + 
Yi xc where Ho is the Hartree Hamiltonian. Eq. 
imposes the locality condition and is equivalent to the 
self-consistency condition of the dynamical mean-field 
theory*^ 



III. CONDUCTION BAND SPECTRA OF 
TRANSITION METALS 

Using the technique described above (below referred 
to as "GW") we calculated the spectral densities N(e) — 
TrlmG(e) for elemental 3d and id metals, where LDA 
has long been the only appropriate approximation. We 
used the fee structure for hep elements (Sc. Ti, Y, Zr, 
Tc, Ru, Co) and the bec structure for Mn. All atomic 
volumes were taken from experiment. Cr was treated 
within the single bec cell, and hence was non-magnetic. 

The iterations were started from the LDA potential, 
and at each iteration the non-local self-energy was mixed 
to this potential with a gradually increasing weight, so 
that in the end only H xc was left. In the final state the 
magnitude of T, xc differs by about 40% from its initial 
value based on the LDA wave functions. The iterational 
procedure is rather stable in all metals except Ni where 
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the magnetic moment is very sensitive to the details of 
the calculation. 

The results are shown in Figs. 1,2 along with the LDA 
densities of states (DOS) and the experimental XPS and 
BIS spectra taken from Refs. Il3lll4l Strictly speaking, 
comparison with experiment requires the calculation of 
the corresponding matrix elements, but we believe that 
in the present context some qualitative conclusions can 
be drawn based on the spectral densities. 

The main differences between the LDA DOS and GW 
spectral density may be summarized as follows. The con- 
duction band widens as all DOS features are pushed away 
from Ef] this outward shift is roughly proportional to 
the distance from Ef- Moreover, all DOS features are 
increasingly smeared out due to the decreasing quasipar- 
ticle lifetime as the distance from Ep increases. Substan- 
tial spectral weight is transferred from the quasiparticle 
states to the incoherent "tail" extending far below Ep. 

The LDA DOS for d-metals is typically too small to 
account for experimentally measured electronic specific 
heat C = jT. As seen from Figs. 1,2, the GW spectral 
density at Ep is generally smaller compared to LDA. 
Although to obtain 7 we have to remove the renormal- 
ization factor Z from ImG, the GW method does not 
improve the overall agreement with 7 measurements. 

The differences between the GW and LDA spectra are 
most notable for early transition metals (treated here in 
the fee structure) where the strongest unoccupied peaks 
given by LDA are almost completely smeared out in GW. 
Notably, the worst agreement with experiment for the 
LDA spectra is observed for the same elements. In par- 
ticular, the position of the Fermi level obtained in LDA 
for hep Sc and Y appears to be off by nearly 1 eVii^ 

The agreement between LDA and GW improves as we 
progress to later transition metals and as the fee structure 
is replaced by the bec one. Already in V and Nb the GW 
and LDA curves are quite similar except for the shift of 
the unoccupied d-peak by 1.5-2 eV to higher energies. 

As in LDA, the 3d metals from Mn to Ni were mag- 
netic in our calculation. For simplicity, we used the bee 
structure for Mn, and also the fee structure for Co (this 
structure is often realized in thin films). As it is seen from 
Fig. 1, the general features of the GW spectral density 
described above are observed in these metals as well. In 
general, the GW description also gives a reasonable ex- 
change splitting and magnetic moment. We obtained the 
moments of 0.9 hb for Mn compared to 1.03 /is in LDA, 



2.3 hb for Fe compared to 2.25 {1b in LDA, and 1.85 [is 
for fee Co compared to 1.62 hb in LDA. From these re- 
sults it is clear that there is no systematic trend for GW 
to give larger or smaller magnetic moments compared 
to LDA. The magnetic moment in Ni is rather sensitive 
to various details of the calculation, and proper conver- 
gence turned out to be problematic. We believe that the 
approach based on Matsubara Green's functions is neces- 
sary to stabilize this problem. Apart from the exchange 
splitting, the shape of the spectrum is quite stable. We 
also note that the magnetic moment is expected to be 
sensitive to the choice of the skeleton graph set. 

As an example of the general trend of band dilatation 
off the Fermi level, the distance from Ep to the upper 
edge of the fully occupied d-band in Cu and Ag is no- 
tably larger in GW compared to LDA, which results in 
a better agreement with experimental XPS spectra. On 
the other hand, we also observe a rather strong upward 
shift of the unoccupied peak in V, Nb, Fe and Mo in ob- 
vious disagreement with the BIS spectra. This shift is 
more pronounced compared to the downward shift of the 
occupied states at a similar distance from Ep. 

The results presented above demonstrate that self- 
consistent GF approach with one-site approximation can 
be used for ground state studies. For 3d and 4d systems 
we have shown that the GW spectral density is generally 
similar to the LDA density of states, while the GF ap- 
proach includes typical Fermi-liquid effects such as finite 
quasiparticle lifetime and self-consistent renormalization. 
The preliminary comparison with experiment is rather 
satisfactory and clearly indicates the problems that need 
improvement beyond our simple GW-OSA technique are 
similar to those in LDA — the exchange splitting in fer- 
romagnets, value of N(Ep), and the unoccupied peak po- 
sition which is too high for some metals. In general, the 
presented technique based on the Luttinger-Ward func- 
tional is a practical alternative to other methods based on 
the density functional theory and can serve as a reliable 
starting point for more sophisticated methods.. 
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FIG. 1: Spectral densities obtained in GW-OSA (red lines), LDA densities of states (black lines), and experimental XPS and 
BIS spectra (blue dotted lines shifted up) for 3d metals. For magnetic elements the calculated curves for majority-spin and 
minority-spin electrons are plotted in solid and dotted lines, respectively. Elements for which the crystal structure was not 
experimental (including Cr) are marked by a star in the upper left corner of the graph. In all graphs the x axis shows energy 
referenced from Ep in eV, and the y axis denotes the spectral density in eV _1 for the calculated curves. 




FIG. 2: Same as Fig. 1 but for 4d metals. 
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